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Abstract 

The Micro-Arcsecond X-ray Imaging Mission (MAXIM), a proposed concept 
for the Structure and Evolution of the Universe (SEU) Black Hole Imager mission, is 
designed to make a ten million-fold improvement in X-ray image clarity of celestial 
objects by providing better than 0.1 micro-arc second imaging. Currently the mission 
architecture comprises 25 spacecraft, 24 as optics modules and one as the detector, which 
will form sparse sub-apertures of a grazing incidence X-ray interferometer covering the 
0.3-10 keV bandpass. This formation must allow for long duration continuous science 
observations and also for reconfiguration that permits re-pointing of the formation. To 
achieve these mission goals, the formation is required to cooperatively point at desired 
targets. Once pointed, the individual elements of the MAXIM formation must remain 
stable, maintaining their relative positions and attitudes below a critical threshold. These 
pointing and formation stability requirements inpact the control and design of the 
formation. 

In this paper, we provide analysis of control efforts that are dependent upon the 
stability and the configuration and dimensions of the MAXIM formation. We emphasize 
the utilization of natural motions in the Lagrangian regions to minimize the control 
efforts and we address continuous control via input feedback linearization (IFL). Results 
provide control cost, configuration options, and capabilities as guidelines for the 
development of this complex mission. 
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I. Introduction 


The MAXIM concept for NASA's Black Hole Imager mission utilizes interferometric techniques 
at the short wavelengths of X-rays to resolve objects as infinitesimal as the Earth at the center of 
the Milky Way, the solar disk in another galaxy, and at greater distances yet, the huge black holes 
that power quasars. While MAXIM's image of a black hole will certainly be one of the new 
century's most spectacular scientific pictures, it will also address leading questions in black hole 
research. In order to achieve these high-precision angular resolution images, very long optical 
baselines are needed. Multiple free-flying spacecraft comprise the sparse aperture providing 
sufficient collecting area as shown in Figure 1. Other components of the interferometric space 
system are a combiner spacecraft in some cases and the detector spacecraft.'’ 2 

Images are generated through interference patterns gathered from the multiple satellites housing 
the optical elements that form the aperture. The interference patterns or fringes are observed only 
if the pathlengths are controlled to great precision (nanometers in some cases). The challenge is to 
control this pathlength in the presence of environmental and spacecraft disturbances driving the 
need for active control systems. The sparse aperture array is configurable since the formation of 
the satellites can change, therefore providing greater sensitivity to a range of angular scales 
allowing the mission concept to meet a variety of target science goals. The science community 
realizes the importance of space-based interferometry and the potential for significant scientific 
breakthroughs. The National Academy of Sciences recommended the development of space- 
based imaging interferometery in its Decadal Report Astronomy and Astrophysics in the New 
Millennium. 

In addition to the Black Hole Imager mission, NASA has several space interferometers planned in 
the future to explore the universe, image stars, detect gravitational waves, and search for life: 
Stellar Imager(SI), Submillimeter Probe of the Evolution of Cosmic Structure (SPECS), Space 
Interferometer Mission (SIM), Terrestrial Planet Finder (TPF), Laser Interferometer Space 
Antenna (LISA), and Planet Imager. Since several of these missions involve formation flying, 
research on formation control of space-based interferometers is extremely important to achieving 
the bold and compelling science waiting to be discovered. 

The concept of formation flight of multiple spacecraft near the libration points of the Sun- 
Earth/Moon (SEM) system offers as many possibilities for space exploration as technical 
challenges, of which the MAXIM mission is a prime example. MAXIM orbits the L 2 Lagrange 
point and necessitates a very stringent line-of-sight (LOS) pointing that affects the formation 
control scheme and the overall configuration. 2,3 That is, the alignment of the telescope array LOS 
(comprised of the separated spacecraft optics and detector) with respect to the target must have 
stability better than the angular resolution during an observation. In the baseline design, the 
individual s/c only require a pointing accuracy of 10 arcsec. We must detect spacecraft 
translations relative to the LOS defined by the target and a fiducial point on one spacecraft. 
Control of the formation at a micron level is required. 

We focus on the dynamics and control of formation flight in a full ephemeris modeling of the 
libration orbit to incorporate all gravitational perturbations. Solar radiation pressure is also 
included. A continuous control technique is applied to control the MAXIM formation to the 
required specifications. The rich natural dynamics near the libration points are explored in depth 
to aid in the development of more complex formations that meet MAXIM formation 
requirements. Our analysis focuses on the amount and duration of the formation control effort 
versus science observation requirements as measured in the formation optics plane. 
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Figure 1 - MAXIM Baseline Mission Concept using Grazing Incident Optics 


A. Modeling assumptions 

MAXIM is comprised of three formation components; the hub spacecraft, ffeeflyer (periscope) 
spacecraft, and the detector spacecraft, which is at a distance of up to 20,000 km from the hub. 
The spacecraft areas and masses are taken from a GSFC Integrated Mission Design Center 
design. The wet mass (bus plus payload) for the Hub is 331 kg, the ffeeflyer spacecraft is 304 kg, 
and the detector spacecraft is 619 kg. The hub and ffeeflyers have the same hexagon shape of 1.3 
x 2 meters while the detector spacecraft is 1.5 x 1.5 meters but has a varying area depending on 
its profile (top view is 1.9 m 2 , side view is 3.3 m 2 ; back/front view is 5.6 m 2 ). 

The science observation duration is 1 e5 secs and is considered a "typical” exposure time required 
to image supermassive black holes. The control effort considers requirements for observations 
with and without thrust accelerations. The controller options for analysis are off during 
observation durations and on in between to realign and maintain the formation. The other option 
is continuously on during observations. The control effort emphasizes continuous control but we 
review a discrete control example as well. 

The targets proposed include the Galactic Center, M87, NGC 4594, Cen A, NGC 4258, NGC 
3227, NGC 5548, NGC 4151, and MCG-6-30-15. While the selection of the targets may affect 
the control effort, a mid point of +45 degree elevation and- 45 degree azimuth as measured in the 
inertial system was chosen for analysis of the maintenance cost. 

The size of the L 2 libration orbit used in the analysis is of a typical mission such as the Wilkinson 
Microwave Anisotropy Probe (WMAP). While the obvious goal is to provide a minimum AV 
budget for the transfer and injection implying a large Lissajous orbit, the effect on the formation 
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control is not as pronounced. Solar radiation pressure is included as the spacecraft are not all the 
same area and mass. 

As seen in Figure 1, the hub and freeflyers form a physical configuration defined here as the 
optics plane which is perpendicular to the detector-hub line of sight (LOS) to a target. The use of 
the optics plane associates a physical configuration (the states of the freeflyers relative to the hub) 
to science requirements derived from a Fourier transform of the image plane, the UV plane. The 
phasing of the mirrors in the ffeeflyer spacecraft obtain a common path length to perform 
interferometry. Maintaining this optics plane physical configuration is our driving requirement for 
the formation control. During one exposure, the distance of an individual ffeeflyer from the hub 
must change by less than 5 microns, so the light contributed by each ffeeflyer to the image moves 
by less than the width of the point spread function. What is important here is that motions of more 
than 5 microns translate into a change in optical path length for that channel of more than 1/1 Oth 
of a wavelength. The MAXIM formation specifies that 25 spacecraft are to be used. Our 
preliminary analysis uses the detector, the hub, and six freeflyers at the maximum nominal radial 
distance of 500 meters from the hub in the optics plane. The current science concept achieves the 
required angular resolution of 0.05 micro-arcseconds with an effective collecting area of 1,000 
cm 2 . This is done by distributing the periscope mirrors among the 25 formation flying spacecraft, 
across a virtual aperture of one kilometer diameter with focal length of 20,000 km and 5 micron 
pixel size. 

The separation between the detector and the hub is also important. For our analysis we used the 
largest distance specified of 20,000 km. For the smaller separation of 200 km, the hub and 
freeflyers would be 'attached', yielding a completely different control effort. A minimal thrust 
level of one pN poses some current technology challenges but may be required to meet the 
assumptions and requirements. 


B. MAXIM Libration Orbit 

The nominal orbit used in our 
analysis places the hub spacecraft on 
a reference that evolves along an L 2 
Halo orbit that has amplitudes of 
700,000 km in Ay and 200,000 km 
in Az. The reference Halo orbit was 
computed with a full Ephemeris 
model and is shown in Figure 2. Six 
freeflyers are evenly spaced, each at 
500-meter distance from hub, on the 
optics plane. The detector spacecraft 
is 20,000 km aft of the hub, along 
the -W direction as seen in Figure 3. 

The target is at +45 deg azimuth and 
+45 deg elevation relative to 
EMEJ2000 Inertial Frame. Solar 
Radiation Pressure (SRP) effects 
based on equal effective areas for 
each spacecraft (2.0106e-6 km2) are Figure 2. MAXIM L2 Halo Orbit with Ay = 700,000 
used. Previous analysis km and Az = 200.000km 

demonstrated that the AY cost for 
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formation maintenance varies over the orbit as a function of the elevation and azimuth angles. 5 
The target angle represents a location that is neither a worst case nor best case for the AV cost to 
maintain the direction. 


C. MAXIM Frame Definition 

The definition of the frame used for our 
control is shown in Figure 4. The X, Y, Z 
axes show the ephemeris inertial axes and the 
local frame is represented by u, v, w. The 
MAXIM hub spacecraft is located at the 
X,Y,Z origin and the angles a, 5 provide the 
alignment toward the target. The target LOS 
is parallel to the w axis. The alignment is from 
the detector to the Hub to the target. 

The local frame components are 
constructed from the following direction 
cosines as 

w = C a C s X + S a C s Y+S s Z 


. Zxw 

U=TZ -1 


u — - 

Zx w| 


A A 

V = wxu 



Figure 3. Orientation Of Hub On Reference Halo 
Orbit and Freeflyer Alignment (not to scale) 


Where the direction cosine matrix (where C=cosine and S = sine) can be constructed as 


'C 


u 


S a C a S s C a C s 
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Figure 4. Maxim Definition for the Inertial (X,Y,Z) and Optics Frames (u,v,w) 
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II. Orbit Control Strategies 


Our investigation takes a global view of the large-scale formation flying problem. Much of the 
research to date considers the continuous control of constellations, clusters, and formations for 
missions near the Earth, or where the influence of other gravitational perturbations can be safely 
ignored. 4 ' 7 Impulsive schemes have also been implemented but these only apply to formations 
that do not require close tracking of a reference solution. 8 Nonlinear methods, combined with 
adaptive control, have also been developed in the two-body problem. 9 ' 10 In contrast, the 
documented research on formation flight in multi-body systems is scarce. Most of these 
investigations focus on the simplified circular restricted three-body problem (CR3BP). 11,12 Some 
authors do consider formation flight relative to Lissajous trajectories developed in the ephemeris 
model but the effectiveness of such a controller is only demonstrated relative to the linear 
dynamics, not the full nonlinear system. That is, the evolution of the controlled formation is not 
numerically integrated in the nonlinear system but rather approximated from the linear dynamics 
of the integrated Lissajous trajectory. 

Previous work by Howell and Barden and others considered naturally occurring formations 
derived from center manifold analysis, as well as a discrete impulsive control approach to 
maintain a prescribed formation plane formed by multiple spacecraft. 11 ' 13 ' 22 One item of interest, 
which pertains to the continuous control approach, is to obtain a rough analytical approximation 
of this center manifold motion and determine how continuous optimal control and exact feedback 
linearization compares, in terms of cost, to the discrete station-keeping approach. 2122 These works 
also demonstrate the efficiency and cost effectiveness of both input feedback linearization (EFL) 
and output feedback linearization (OFL) methods for formation control in the CR3BP. The IFL 
controller is designed to force the error dynamics of each state variable to follow a critically 
damped response. The OFL controller, on the other hand, was only initially designed to force the 
radial separation between each spacecraft to track some specific value. Hence, no relative 
orientation requirements were imposed on the formation. But, the OFL controller can also 
constrain the orientation rate and that may be beneficial in certain applications. A linear quadratic 
regulator (LQR), derived from optimal control theory, yields essentially an identical error 
response and control acceleration history as the input feedback linearization approach. However, 
the IFL controller is computationally much less intensive and, by comparison, conceptually 
simple. We address the properties of the IFL controller in defining the MAXIM formation 
control. 

Our analysis considers strategies to maintain a planar formation of the spacecraft in an orbit about 
the Sun-Earth L 2 point. That is, controlling the deviations of each spacecraft relative to the initial 
formation plane. For a comparison, a discrete stationkeeping control approach is devised to force 
the orientation of the formation plane to remain fixed inertially. 


A. Discrete or IFL Control 

Even if the onboard propulsion system can deliver an accurate continuous thrust level, it is likely 
that the science requirements of the mission may stipulate some thruster down-time. This leads 
then to a new consideration about the degree of accuracy to which the formation can be 
maintained via discrete impulses applied at the beginning of each N segments. A simple targeter 
applying discrete impulses can maintain orientation. Alternatively, an LQR controller, based on a 
discrete time system, yields the optimal magnitude of each differential control impulse at the 
specified time interval. However, in either case, the nominal control input that must be added is 


6 

American Institute of Aeronautics and Astronautics 



still assumed to be available as needed. If a continuous nominal control effort is not available, 
how often must impulsive maneuvers be implemented to maintain the nominal configuration? 

Of the two options mentioned here, the most straightforward approach is a simple strategy that 
targets the state vector at the end of each specified interval. Consider a case of a 100 km 
formation of two spacecraft, constrained to remain aligned with the inertial y-axis at all times. In 
the absence of any external perturbations, a continuous control approach via IFL requires 0.3348 
m/sec of total AV over a period of 180 days. How much will the formation diverge if the control 
effort is discretized over a period of hours, or even days? Consider the general form of the 
solution to the linear system. 


Sj\ 


Sr 0 


'A 

B~ 

Sr, 1 



— 

C 

D 

_<Sv 0 '+Av 0 _ 


where <t>(r, ,/ 0 ) denotes the state transition matrix associated with the nominal Lissajous orbit 

along which the MAXIM detector spacecraft is assumed to evolve. The symbol S denotes a 
perturbation relative to the nominal Lissajous trajectory and Av 0 represents an impulsive 
maneuver applied at the beginning of a segment. The superscripts identify velocity perturbations 
before and after a maneuver. Controlling the position of the formation spacecraft relative to the 
detector or hub to a constant vector, as observed in the inertial frame, is equivalent to targeting a 
particular constant perturbation 5r x relative to the inertial frame. An impulsive maneuver of the 
form 

A v 0 = B-'(Sf[-ASr 0 )-5v- 

will accomplish die goal in the linear system. A precise implementation of this scheme in the 
nonlinear system requires that the maneuvers be differentially corrected over each segment. 
However, the simple expression is found to accomplish the goal, provided the length of each 
segment is sufficiently small. 

A discrete optimal control strategy was also developed and applied to this problem and was 
formulated to implement impulsive maneuvers at the beginning of each segment. The cost 
function was defined as the sum of two terms. The first term was a weighted (W) quadratic 
function of the state errors at the end state; the second term is a weighted (Q) quadratic function 
of the state deviations along the path. Of course, the simple corrector only targets the states at the 
end of each segment so the deviations of the six states along the path will be greatest with a 
discrete set of impulses determined via such a targeting scheme. Along the path, the optimal 
controller was formulated to weight position deviations to be twice as costly as velocity state 
errors. Then, the end-state weighting can be varied. 

Not surprisingly, the total AV computed via either discrete control approach converges on the 
cost computed via the continuous control approach as the interval between maneuvers becomes 
shorter. The state corrector only targets position, not velocity. Hence, the rate of divergence is 
increased because the velocity error at the terminal point of each segment is not zero. This 
position divergence appears for a MAXIM example in Figure 5 for a sampling of three distinct 
discrete controllers. The dashed red line in the figure is a simple targeter. At a given time, an 
impulse targets the required position at a time in the future, in this case 6 days. The plot shows 
the position deviations between maneuvers. Note that if we cut the time interval in half, we will 
cut the deviations in half. The blue is a discrete optimal controller that minimized the state 
deviations without end point weighting. (Note that the largest amplitude points along the blue 
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Nominal Position Vector = 500-m Along Inertial Y-Axis 


curve correspond to the impulses.) This optimal controller does cut the state error between 
maneuvers. The green line is generated as a discrete optimal controller but with an increase in the 
weighting on the end state (Q m = le' 2 
Q). The deviations naturally increased. 

If the end-state weighting increases a 
sufficient amount, the deviations are 
the same as the simple targeter. 

Actually there is a fourth curve in the 
figure. A magenta curve is plotted for 
the optimal controller with a large 
weighting on the end states. It is, in 
fact, seemingly identical in the figure 
to the red curve produced with the 
simple targeter. In any case, once the 
spacecraft reaches the target position it 
will move away from it because the 
relative velocity is not zero. Note that 
the deviation is in meters and implies 
that to meet the MAXIM requirements, 
the discrete control would need to run 
near continuously. An interval of only 
a few seconds is available before the constraints are violated. This leads us to use a true 
continuous controller, that is, the input feedback linearization method. 



Time (days) 

Figure 5. MAXIM Position Deviation with Discrete 
Control with Maneuvers Everv Six Davs 


In general, the control goal should seek to minimize the integral of the error while only resorting 
to discretized maneuvers, which ensures that the overall error is minimized everywhere along the 
segment rather than at the endpoints alone. This minimization is one of the issues to be further 
addressed in this continuing investigation. Another aspect of this study considers the robustness 
of these control strategies in the presence of modeling uncertainties. 

B. Control Methodology 

It was first thought that the control methodology might be based on the use of a paraboloid as a 
control surface for the freeflyers as shown in Figure 6. 19 Here the hub is located at the center of 
the paraboloid and the parabola shape constrains the freeflyer spacecraft. This geometric shape 
allows for a controller that can accommodate freeflyers at various radii and has been used 
successfully to shift freeflyers from one radius to another. When the problem is formulated in 

w (LOS) 

Zenith 

x= a pJ^JK cosv 

y = a pyl u p /h p sinv 

z = u p +q 

— HD ~ /v ~ /\ ~ /\ 

r E ' =xu + yv + zw 


Figure 6. Possible Parabolic Configuration Concept 
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terms of the paraboloid variables 
(q,v,u ), however, the controller 

is not designed for cases when the 
detector or any of the freeflyers 
are along the line of sight. The 
difficulty is a singularity in this 
geometry. In this configuration, 
the freeflyers (deputy spacecraft) 
are placed on the dashed green 
line and evenly spaced. Thus, 
alternatively, the control problem 
is formulated in Cartesian 
coordinates that are related to the 
paraboloid variables as needed. 



Figure 7 Nominal Motion and Determination of Vehicle 
Position Relative to UVW-Frame 


This formulation can be modified to the configuration shown in Figure 7 where the angle v is the 
parameter of importance for science observations. In Figure 7, the nominal behavior would 
correspond to e = 0. The angle v is measured in a counterclockwise direction (similar to azimuth) 
Each freeflyer has its own d x and the kinematics are written as 

r m =rd x 

u r HD ' -rd x + rvC £ d 2 + rsd 3 
d x — C\C v u + C c S v v+ S £ w 

We define a conservative displacement tolerance in v and of 5-microns. Each freeflyer is at a 

given nominal orientation with respect to the hub with a value of v and e. To relate the Cartesian 
coordinates to spherical, 

x = rC v C £ x = rC v C £ - rvS v C £ - reC v S £ 

y = rS » C s y = rs v c £ + rvC v C £ - reS v S £ 

z = r S e z = rS £ + reC £ 

The nominal motion is in the spherical coordinates while the control effort is formulated in the 
inertial focal frame. Using initial conditions of a 1 km diameter plane for the freeflyer spacecraft, 
the formation is defined in terms of the spherical coordinates as follows. For this example, it is 
easily translated into the Cartesian coordinates. Note again that these are components in and 


normal to the optics plane. 




r = 20,000 km 

x =0 

| x=rC, 

r = 0 

v = 0 

VII; 

II 

o 

Y y= 0 

v> = 0 

£ = —7lj2 

z= 0 

1 z* =rS £ 

£ = 0 

r =1 km 

r =0 -'j 

II 

* 

x=0 

v* = n/r/3;(n = 0...5) 

v* = 0 

y = rS . 

V 

O 

II 

• 

£*=0 

£ =0 J 

9 z’ =0 

o 

II 


American Institute of Aeronautics and Astronautics 








C. Controller Development 


20-22 


The controlled equations of motion (EOM) in the ephemeris frame (as related to an Earth 
centered frame) are expressed where P 2 is the Earth. Here “D” can be representative of either the 
deputy (ffeeflyer) or the detector in terms of using the equation for control. The superscript on 
the left denotes the frame used for formulation of the differential equations. 




(- 3 A ijrPiD, 


-/(*> 
'f, h “ =J(r, r '“ 


)+«! D ' , (0 


Since the hub is assumed to move in a Lissajous orbit computed in a full ephemeris model, no 
control is included for its evolution along this path. The controller is selected to the type of 
response that is required, and in most cases this is a critically converging system. The above 
terms are subtracted to yield the EOM for the relative motion in both the inertial and the u,v,w 
frame by using the direction cosine matrix for rotation into the focal plane. The controller in the 
local frame is constructed as follows, 


% m ‘ ={"C'} Aj, + {“C'} sf■> (<) = { 1 "C'} AT, + S«> (() 


Here the dynamics of the system can be annihilated by the first term in the IFL controller and the 
remaining terms show the critical response chosen for the IFL controller. 

= -r)-<o n 1 (r u m '-7) 

-► {'cT”^=a/,+*,«>(,) 

Once the control is determined in the optics frame to generate a critical response, it is translated 
back into the inertial frame for use by the control algorithm, 

M / (Z) ' ) (f) = { / C C, }M (D ' ) (f) 

It is assumed that this control is available continuously." 


m. RESULTS 

Figure 8 shows the orientation of a sample MAXIM formation. The left view is the optics plane 
view and the arbitrary positions of the six freeflyers relative to the hub. On the right is the inertial 
frame view. The cost over one revolution of a halo or Lissajous for each freeflyer and the detector 
is slightly different for each but generally similar. In these cases the solar radiation pressure is 
turned on. An effective area of 1.6m 2 is assumed. 
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UV-Plane View Inertial Frame View (a=45°,5=45°) 



u (km) 


Figure 8. MAXIM formation design for IFL control 


A. Maintenance Maneuvers 

Maintenance maneuvers must be performed during the mission to preserve the strict separation 
and orientation requirements. Using IFL, a control effort for 180 days (one libration orbit 
revolution) was computed and shown in Figure 9 are the thrust profiles for both the detector and 
the freeflyer spacecraft. The effort is on the order of several mN for the detector and tenths of pN 

Detector Thrust Profile 



Free Flyer Thrust Profiles 
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Control Input (p N) Control Input (mN) 


for the freeflyers. If the mass ratios are changed, doubling the mass, the thrust profile changes 
accordingly requiring twice as much thrust as shown in Figure 10. 



B. Maintenance and Recovery 

We considered the effect on turning the thrusters off during a MAXIM observation, which is 
100,000 seconds long in duration. This was done to analyze any requirement of eliminating 
accelerations from the small but continuous maintenance. This also allowed us to investigate the 
drift in the formation over a small time period and to establish the thrust profile necessary to 
recover the original formation alignment. Figure 11 shows the increase in radial errors of the 
detector and the freeflyers with respect to the hub over the observation duration. Once the 
100,000 seconds is reached, the controller is activated again to reset the freeflyer and detector for 
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the nominal configuration, without any error. Note that the errors are not linear and ramp up in 
the typical exponential fashion for a libration orbit. The peak radial error of 15 km was observed 
for the detector at a distance of 20,000 km from the hub. The freeflyer radial errors were smaller, 
ranging from 300 mm to 550 mm. Figures 12 and 13 show the enrors in the optics frame (u,v,w) 
components. Also apparent is the in-plane (u and v components) and out-of-plane deviation. 
Figure 14 shows these enrors as pointing errors in arc-seconds. The formation requirements for 
the freeflyers includes a radial error in the optics plane of less than five microns, clearly less than 
1mm, and an azimuthal angle error (v-v*) less than 5 micro-arcseconds. One can easily see that 
the drift results in the science requirements being violated. It is also apparent that the radial error 
exceeds the requirement very quickly and shortening the duration is not a feasible solution. 



Taneidayt) 

Figure 12, Freeflyer Focal Plane Error Figure 13 Detector Focal Plane Error 

(u,v,w) during No Control 
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Figure 15. Three Day Span with Maintenance, No Control, and Recovery 

To assess the thrust level and compare on-off thrust magnitude, Figure 15 shows the thrust profile 
for the detector and freeflyers during a three day span that includes 1 day of maintenance, 
100,000 seconds of no-control due to science observations, and the recovery over less than !4 day. 
The recovery thrust levels reached one N for the detector (top plot) while the freeflyers (bottom 
plot) required less than 10 pN. Figure 16 shows the thrust profile for the maintenance (segment 
A) of the 3-day simulation. The detector control required 3e-3 N while the freeflyers required less 
than 0.05 pN. During maintenance the scale of the plot implies that the thrust level is constant. In 
reality, even though the control offsets the perturbations to produce zero error, the thrust level 
varies because the forces are not constant. This will be more apparent when errors are added to 
the model. During observation, i.e., no control (segment B) appears as expected and the thrust 
profile is zero. The recovery thrust profile corresponding to segment C appears in Figure 17. 



Time (days) Time (days) 

Figure 16. Maintenance Thrust Profile Figure 17. Recovery Thrust Profile 

for Segment A for Segment C 
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C. Reorientation 


Another important aspect of the mission is reorientation. In this preliminary investigation, it was 
also of interest to determine the ability to reorient and the relative cost for a simple scenario. This 
simple scenario also avoids any consideration of collisions at this time. Reorientation of the 
formation to view additional targets was accomplished by performing a 90-degree rotation about 
the inertial z-axis as shown in Figure 18. Initially, the target was assumed to be along the inertial 



Figure 18. Formation Reorientation, 90 Degree Change in X-Y plane 


position x-direction. This x-axis pointing was then re-aligned to the y-axis. The elevation angle 
was held to zero degrees. As seen in Figure 18, the initial to final orientation shows the switch in 
the direction of the axes with respect to the freeflyer spacecraft. The initial nominal angle 
information for v and £ for this sequence is shown in the bottom left. 

The thrust levels for the detector and the freeflyer are plotted in Figure 19 over the duration of the 
reorientation. The thrust profile for the reorientation was allowed a duration of 7 days to 
minimize the effort of the controller for the detector spacecraft. Even with this amount of time 
the detector thrust level reached a maximum of 1.5 N and took approximately 8 days to 
completely recover and settle into the proper configuration. Of course, in this example, the 
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detector is 20,000 km away. The thrust profile for the ffeeflyers was much lower at values less 
than 25pN. Note the secondary increase in the thrust magnitude. This is an artifact of the control 
as the thrust compensates for the velocities that must also be attained for the proper final state 
vector of each spacecraft. In Figure 20, the reconfiguration is represented in ephemeris inertial 
coordinates. Note the switch in the x and y components that show the repositioning of the 
spacecraft for the reorientation. 






Figure 20. Thrust Profile for Reorientation of 90 degrees 
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IV. Summary 


Two approaches, discrete and continuous, were investigated for the control of the Maxim 
formation. It was found that control of the MAXIM formation can be accomplished using either 
discrete control with a simple or an optimal control strategy or by Input Feedback Linearization 
(IFL) control. To maintain the 5-micron tolerance on the freeflyer positions, the discrete control 
needs to operate at a time interval that approaches a continuous effort. An IFL was used which 
combines the effect of annihilating the environmental dynamics while adding a specific user- 
defined critically damped response to achieve the maintenance accuracy or to perform a recovery 
within a reasonable time. 

The total effort of maintaining the nominal formation in the halo orbit requires a detector thrust 
level that ranges from 4 mN to 7 mN and freeflyer thrust levels of 0.1 pN to 0.3 pN. If allowed 
to drift, the science requirement of maintaining a five-micron position of the freeflyer in the 
optics plane is almost immediately violated, indicating the need for continuous thrusting. The 
thrust level for formation recovery after a reorientation ranges from between 0 to 1 N for the 
detector to less than a maximum of 15 pN for the freeflyers. These efforts do not include 
navigation or maneuver errors or navigation measurement updates. 

V. Conclusions 

It was found that a continuous control effort is required to maintain the MAXIM formation 
requirements. The minimum thrust levels are at pN levels, which can be achieved by current 
propulsion technology. The challenge is whether a propulsion system can be implemented with 
control at the pN level and with the required power to maintain the formation if no allowable drift 
is permitted. 
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